Yang- Yang thermodynamics on an atom chip 
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We investigate the behavior of a weakly interacting nearly one-dimensional (ID) trapped Bose 
gas at finite temperature. We perform in situ measurements of spatial density profiles and show 
that they are very well described by a model based on exact solutions obtained using the Yang- Yang 
thermodynamic formalism, in a regime where other, approximate theoretical approaches fail. We 
use Bose-gas focusing [Shvarchuck et al., Phys. Rev. Lett. 89, 270404 (2002)] to probe the axial 
momentum distribution of the gas, and find good agreement with the in situ results. 
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Reducing the dimensionality in a quantum system can 
have dramatic consequences. For example, the ID Bose 
gas with repulsive delta-function interaction exhibits a 
surprisingly rich variety of physical regimes that is not 
present in 2D or 3D [3,13]. This ID Bose gas model is of 
particular interest because exact solutions for the many- 
body eigenstates can be obtained using a Bethe ansatz 
0| . Furthermore, the finite-temperature equilibrium can 
be studied using the Yang- Yang thermodynamic formal- 
ism 0, [H, Q , a method also known as the thermodynamic 
Bethe ansatz. This formalism is the unifying frame- 
work for the thermodynamics of a wide range of exactly 
solvable models. It yields solutions to a number of im- 
portant interacting many-body quantum systems and as 
such provides critical benchmarks to condensed-matter 
physics and field theory [f|. The specific case of the ID 
Bose gas as originally solved by Yang and Yang Q is of 
particular interest because it is the simplest example of 
the formalism. The experimental achievement of ultra- 
cold atomic Bose gases in the ID regime Q has attracted 
renewed attention to the ID Bose gas problem Q and is 
now providing previously unattainable opportunities to 
test the Yang- Yang thermodynamics. 

In this paper, we present the first direct comparison 
between experiments and theory based on the Yang- Yang 
exact solutions. The comparison is done in the weakly 
interacting regime and covers a wide parameter range 
where conventional models fail to quantitatively describe 
in situ measured spatial density profiles. Furthermore, 
we use Bose-gas focusing [§] to probe the equilibrium 
momentum distribution of the ID gas, which is difficult 
to obtain through other means. 

For a uniform ID Bose gas, the key parameter is the 
dimensionless interaction strength 7 = mg/h 2 n, where m 
is the mass of the particles, n is the ID density, and g is 
the ID coupling constant. At low densities or large cou- 
pling strength such that 7> 1, the gas is in the strongly 
interacting or Tonks-Girardeau regime [r|. The oppo- 
site limit 7 <C 1 corresponds to the weakly interacting 



gas. Here, for temperatures below the degeneracy tem- 
perature Td = h 2 n 2 /2mkB, one distinguishes two regimes 
(i) For sufficiently low temperatures, T <C \f^Tdi 
the equilibrium state is a quasi-condensate with sup- 
pressed density fluctuations. The system can be treated 
by the mean-field approach and by the Bogoliubov the- 
ory of excitations. The ID character manifests itself 
through long-wavelength phase fluctuations resulting in 
a finite phase coherence length 1$ which greatly exceeds 
the mean- field correlation length l c . (ii) The tempera- 
ture interval ^Td <C T <C Td corresponds to the quan- 

where both the density and 
the condition l c <C required 



turn decoherent regime [11 1 
the phase fluctuate. Here 
for the existence of a quasi-condensate is no longer sat- 
isfied and the system can be treated as a degenerate 
ideal gas combined with perturbation theory in g. At 
temperatures near the crossover to the quasi-condensate, 
T ~ ^/jTd, neither of the approximate approaches work 
and one has to rely on the numerical solution to the exact 
Yang- Yang equations, as we show in the present work. 

Experiments on ID gases arc usually carried out in 
harmonic traps with strong transverse confinement and 
weak confinement along the axis, u>± ^> wb. A trapped 
gas is in the ID regime if both temperature and chemical 
potential are small with respect to the radial excitation 
energy, fc^T, \i -C hu>±. It was recognized early on that 
the physics of the degenerate part of the trapped cloud 
is already effectively ID if the weaker condition [i < hui± 
is satisfied Q. For a gas in the ID regime, the effective 
ID coupling can be expressed through the 3D scattering; 
a if a < (ft/row jj 1/2 



2huj 1 



12] 



length a as g 

Various physical regimes of a harmonically trapped ID 
gas have been discussed in Refs. 0, [l3 , 14 , T^j • The above 
classification of the regimes for the uniform gas can be 
applied locally to the trapped gas if the conditions for the 
local density approximation (LDA) are met 0, Q, EH • 
The trapped cloud is then characterized by a global tem- 
perature and a spatially varying density, n — > n(x). 
Thus, a trapped ID gas which is in the quasi-condensate 
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regime in its centre [T <C \Z"f(0)Td(0)] will cross over 
towards the wings - through the quantum-decohercnt 
regime - to the nondegenerate regime [Td(x) <C T] as 
the local density n(x) decreases. 

We experimentally investigate the behavior of a weakly 
interacting trapped Bose gas [7(0) ~ 10~ 2 ] in the regime 
where /i < hu>± and ksT ~ huj_. Similar measurements 
to our in situ data were previously performed at higher 
temperatures and higher linear densities 3, 17|; the ob- 
served density profiles were found to be in disagreement 
with approximate theories. Our approach here is differ- 
ent in that we present a full description of the density pro- 
files through a consistent model based on the solutions 
to the exact Yang- Yang equations [1, EH and that we 
compare real-space and momentum-space distributions 
for the same experimental parameters. 

For (nearly) ID clouds, it is difficult to probe the ax- 
ial momentum distribution. The conventional time-of- 
flight method does not work, mainly because the cloud 
hardly expands axially beyond its initial length. The 
analysis is further complicated by strong density fluc- 
tuations that develop in time-of-ffight from the initial 
phase fluctuations [18j. The axial momentum distribu- 
tion of phase-fluctuating 3D condensates has been mea- 
sured using Bragg spectroscopy , but this technique is 
extremely sensitive to vibrations and requires averaging 
over many realizations. 

We gain experimental access to the axial momentum 
distribution using Bose gas focusing. This is a tech- 
nique introduced for non-equilibrium 3D condensates Q ; 
here, we extend it to study the equilibrium properties of 
(nearly) ID clouds. A detailed analysis of focusing in this 
regime will be presented elsewhere. In brief, we apply a 
short, strong axial harmonic potential yielding a kick to 
the atoms proportional to their distance from the trap 
center (analogous to the action of a lens in optics), fol- 
lowed by free propagation. As a result the atoms come 
to a focus, at which time the axial density distribution 
reflects the axial momentum distribution before focusing. 
Initialphase fluctuations result in a finite width of the 
cloud [9j . Since the focusing brings all atoms together ax- 
ially, the signal level is high, even for a single realization. 
As we will show, averaging over a few shots is sufficient 
to obtain high signal-to-noise ratio. 

The core of our setup is a magnetic microtrap consist- 
ing of three layers of current carrying wires. The sur- 
face layer is formed by a silicon substrate coated with a 
1.8-^m-thick vapor-deposited patterned gold layer. This 
atom chip faces in the direction of gravity (+z). On 
the chip we use a Z-shaped wire, with a 3 mm long and 
125 /Ltm wide central section along x. The remaining two 
layers are behind the silicon substrate and each contain 
three parallel copper wires (diameter 300 /jm), in the x 
and y-direction, centered at z = —0.5 mm and —0.8 mm 
from the chip surface, respectively. The wires in the y 
direction are used to vary the confinement along x. 
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FIG. 1: Linear atomic density from absorption images ob- 
tained in situ (a)-(d) and in focus (e)-(h) by lowering (from 
top to bottom as indicated) the final RF evaporation fre- 
quency. In situ: solid lines are fits using Yang- Yang thermo- 
dynamic equations (see text). The values of fi and T resulting 
from the fits are shown in the figure. Dotted line: ideal Bose 
gas profile showing divergence for n(x) = 0. Dashed line in 
(d): quasi-condensate profile with the same peak density as 
the experimental data. In focus: solid lines are the sum of two 
independent Gaussian fits - one to the wings (dotted lines) 
and one to the central part of the density profile. 



We trap 2 x 10 7 87 Rb atoms in the F = 2,m F = 2 
state in a tight magnetic trap near the chip surface, and 
perform forced evaporative cooling by applying a radio 
frequency (RF) field. The frequency lurf is ramped down 
from 27 MHz to 1.7 MHz relatively quickly (in 180 ms) 
to purposely reduce the atom number, in order to pene- 
trate deeper into the ID regime. Before reaching degen- 
eracy we relax the axial confinement to a final trap with 
lu±/2it = 3280 Hz, w»/2k = 8.5 Hz, and a bottom cor- 
responding to ojuf/Zk = 1.518(2) MHz. The current in 
the Z-wire is set at 2.25 A, and the distance of the cloud 
to the chip surface is 90 /im. In this trap we perform a 
slower ramp (450 ms) to the final RF frequency An addi- 
tional 300 ms of plain evaporation allows the damping of 
residual quadrupole collective oscillations in the cloud to 
the point where these oscillations are no longer visible. 
We probe the gas using standard absorption imaging, 
with a measured optical resolution limit of 4 /im. The 
CCD camera square pixel size is 2.15 in the object 
plane. Our axial trapping potential was characterized 
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using both the measurement of in situ density profiles at 
high T and the dipole mode oscillation frequency in the 
trap center. The curvature in the center corresponds to 
a frequency of 8.5 Hz and gradually decreases to ~ 6 Hz 
in the wings of the cloud [20j • 

In Fig. DJa)-(d) we show the linear density of atomic 
clouds in the magnetic trap for different final RF fre- 
quencies. These data were obtained by in situ absorption 
imaging and integrating the atom number along z. The 
absolute atom number was calibrated using time-of-flight 
data. Each curve is an average of ~ 18 images taken un- 
der identical circumstances. Since all of our data was 
taken for ji < Hlo± [ill ], we expect that the interactions 
will significantly affect only the distribution in the radial 
ground state, while the population in the radially excited 
states can to a good approximation be described by the 
ideal-gas distribution. This leads to the following model 
that was used to analyze the in situ data. 

We start from the solution to the Yang- Yang integral 
equations for a finite-temperature uniform ID Bose gas 
at thermal equilibrium Q ■ This yields numerical results 
for both the equation of state uyy(^,T) and the local 
pair correlation 11 1. The LDA is then used to account 



for the axial potential via a varying chemical potential 
fi(x) = fi — V{x). This approach is expected to be valid 
as long as the axial potential is smooth on the scale 

Since our 
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of the relevant correlation lengths 
temperature is on the order of the radial level splitting, 
hiu^/ks = 158 nK, the fraction of the gas which occu- 
pies radially excited states can not be neglected. We ac- 
count for this fraction by summing over radially excited 
states (radial quantum number j > 1, degeneracy j + 1) 
and treating each radial state as an independent ideal 
ID Bose gas in thermal equilibrium with t he g as in the 
radial ground state, Hj(x) = fJ,(x) — jhio± 22|. Within 
this model, the linear density is given by 

For the radially excited states, we use the result 
of the LDA for the ID ideal gas, n e (fj,j,T) = 
gi/2(exp(iJ 1 j /ksT))/ At where gi/2 is a Bose function 
and A T = (2ttH 2 /mknT) 1 / 2 is the thermal de Broglic 
wavelength [l|, [l5j. Note that as long as fx < fio>j_ we 
have (itj < which is necessary to avoid divergence of 
g±/2- In this model, the radially excited states act as 
a bath for particle and energy exchange with the radial 
ground state. The resulting fits are shown as solid lines 
in Fig. QJa)-(d) and describe our data very well. The 
fitted values of T and \i are displayed in Fig. [2] 

We now turn to the in focus measurements which give 
access to the axial momentum distribution of the gas. 
The focusing pulse is created by ramping up the axial 
trapping frequency from 8.5 Hz to 20 Hz in 0.8 ms, main- 
taining this for 3.8 ms, and ramping back to 8.5 Hz in 
0.8 ms, followed by a sudden switch-off of the magnetic 




FIG. 2: Characterization of the atomic clouds as a function of 
the final RF frequency, as determined from fits of the Yang- 
Yang model to the in situ data and Gaussian fits to the in 
focus data, (a) Temperature from the in situ data (■) and 
from the radial (□) and axial (o) size of the broad Gaussian 
fit to the in focus data. The dash-dotted line is to guide 
the eye and indicates a ratio of 11 of the trap depth and 
the cloud temperature; dashed line corresponds to hu>±/kB- 
(b) Chemical potential from the Yang- Yang fit; dashed line 
indicates huJx/kB- (c) Atom number from the in focus data: 
wide distribution (o) and central peak (A); from the Yang- 
Yang model fit to the in situ data: atoms in the radial ground 
state (■), in radially excited states (•), and atoms in the 
radial ground state experiencing fi(x) > (A). 



trap. During the focusing pulse the cloud length reduces 
by less than 20%. After switching off the magnetic trap, 
the cloud expands in the radial direction on a timescale of 
l/wj_, so that the interactions vanish rapidly compared 
to the relevant axial timescale and the subsequent ax- 
ial contraction can be treated as free propagation. After 
13 ms of free propagation the cloud comes to a focus. 

In Fig.[]Je)-(h) we show the axial density distribution 
obtained in the focus, averaged over typically 10 shots, 
for final RF frequencies similar to the in situ data in 
Fig. [1] (a)-(d). Here, in contrast to the in situ results, 
one can clearly distinguish a narrow peak from a broad 
pedestal for RF values below 1.56 MHz [Figs.QIg)-(h)]. 
The Yang- Yang solution does not yield the momentum 
distribution and thus it can not be used to fit to the in 
focus data. Instead, to quantify the observation of the 
bimodal structure we first fit a 2D Gaussian to the wings 
of the atomic density distribution. In a second step we 
fit a narrow Gaussian to the residual peak in the cen- 
ter. The fitted curves are shown after integration in the 
z direction in Fig. QJe)-(h), and describe the observed 
in focus distributions well. Fig. [He) shows the result- 



4 



ing atom numbers in the wide and narrow part of the 
momentum distribution; we also plot the atom numbers 
from the Yang- Yang model in the radial ground state, 
in the radially excited states, and atoms in the radial 
ground state experiencing y,(x) > 0. Comparing the in 
situ and the in focus data, we conclude that: (i) the mo- 
mentum distribution becomes bimodal around the point 
where the global chemical potential fi crosses zero and 
becomes positive; and (ii) the narrow part of the momen- 
tum distribution is dominated by the atoms in the radial 
ground state (described by nyy), while the wide part is 
dominated by atoms in the radially excited states. The 
shape of the narrow peak thus constitutes a measurement 
of the momentum distribution of the Yang- Yang gas, for 
which currently no theoretical comparison is available. 

A further comparison between the in focus and in situ 
results can be made as follows. Estimates for the tem- 
perature can be obtained from the Gaussian fit to the 
wide part of the in focus data, by assuming that the 
tails (where degeneracy is negligible) are well described 
by Boltzmann statistics. The resulting temperatures are 
shown in Fig. [2Ja). The agreement with the temperature 
extracted from the in situ data is reasonable. We at- 
tribute the remaining discrepancy to the approximations 
implicit in the above interpretation of the Gaussian fit 
results, which neglects the discrete radial level structure 
and the contribution of the radial ground state to the 
wide part of the axial momentum distribution. 

The failure of the ideal-gas and quasi-condensate de- 
scriptions is illustrated in Fig. QJc, d). The key point 
here is the following. The Yang- Yang thermodynamic 
equations yield a smooth equation of state nyy(/i,T), 
including the region around fi(x) = 0. This deviates dra- 
matically from both the ideal-gas description (diverging 
density as fi — > from below) and the quasi-condensate 
description (vanishing density as /i — > from above). 
The region in fi(x) (and consequently in nyy(i)) where 
this discrepancy is significant [la ] is particularly large for 
our parameters, and the Yang- Yang thermodynamic so- 
lutions are essential for a proper description of the data. 

In conclusion, we have found excellent agreement be- 
tween in situ measurements of the spatial linear density 
of a nearly ID trapped Bose gas and a model based on 
the exact Yang- Yang solutions. We have measured the 
corresponding momentum distribution, and expect this 
to stimulate further theoretical work. More generally, 
our results establish a new and strong link between ex- 
periments on low-dimensional quantum gases and exact 
theory for interacting quantum many-body systems at 
finite temperatures, both fields of significant current in- 
terest. Additionally, our findings should be relevant to 
related experimental systems, such as guided-wave atom 
lasers 23 1 and atom-chip based interferometers ■ 
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